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SUMMARY 


Transonic  flutter  of  a  NACA64A00G  airfoil  undergoing  plunging 
and  pitching  oscillations  is  studied  using  the  U-g  and  p-k  methods.  The 
aerodynamic  coefficients  are  calculated  using  an  improved  version  of  an 
ONERA  unsteady  transonic  aerodynamics  code  which  include  the  second 
time  derivative  term  of  the  velocity  potential  in  the  governing  equation. 
Comparisons  with  LTRAN2-NLR  show  good  agreement  up  to  and  in  some 
cases  exceeding  kc  =  0.4,  except  for  the  pitching  moment  curves  at  the 
transonic  dip  Mach  number  of  0.85.  All  flutter  results  are  presented  for 
M  =  0.85.  The  p-k  method  gives  flutter  speeds  identical  to  those  from  the 
U-g  method.  Subcritical  damping  ratios  using  the  U-g  method  with  Frueh’s 
and  Miller's  damping  formula  are  quite  close  to  those  obtained  from  the  p-k 
method,  especially  for  large  values  of  the  airfoil-air  mass  ratio.  Response  of 
the  airfoil  to  externally  applied  forces  and  moments  is  studied  using  the  p-k 
method  and  a  viscous  damping  model  for  coupled  motions. 


RESUME 


II  s’agit  de  l’etude  dcs  vibrations  aeroelastiques  en  regime 
transsonique  d’un  profil  NACA64A006  soumis  a  des  oscillations  dans  le 
plan  vertical  au  moyen  des  methodes  U-g  et  p-k.  Les  coerficients  aerodyna- 
miques  sont  calcules  a  partir  d’une  version  amelioree  du  code  des  vitesses 
aerodynamiques  instationnaires  en  regime  transsonique  de  l’ONERA, 
laquelle  a  recours  a  la  seconde  derivatif  du  temps  du  potentiel  de  la  vitesse 
dans  l’equation  principale.  Des  comparaisons  avec  LTRAN2-NLR  montrent 
que  les  valeurs  obtenues  sont  conformes  jusqu’a,  et  quelquefois  surpassant, 
kc  =  0.4,  sauf  dans  le  cas  des  couibes  de  moment  de  tangage  pour  un  nombre 
de  Mach  creux  transsonic  de  0.85.  La  methode  p-k  foumit  des  vitesses  de 
vibrations  aeroelastiques  identiques  a  celles  de  la  methode  U-g.  Des  rapports 
d’amortissement  moins  critiques  employes  dans  le  cadre  de  la  methode  U-g 
avec  les  formules  d’amortissement  de  Frueh  et  de  Miller  sont  tres  proches 
de  ceux  ohtenus  avec  la  methode  p-k,  surtout  pour  les  grandes  valeurs  de 
rapport  de  masse  de  profile-d’air.  La  reponse  du  profil  a  des  forces  et  a  des 
moments  exterieurs  est  etudiee  par  la  methode  p-k  et  un  modele  d’amortis¬ 
sement  visqueux  pour  les  mouvements  couples. 


CONTENTS 


Page 

SUMMARY . (iii) 

SYMBOLS . (vii) 

1.0  INTRODUCTION . 1 

2.0  FLUTTER  ANALYSIS  OF  TWO-DEGREE-OF-FREEDOM  AIRFOIL  MOTION .  2 

2.1  The  U-g  Method .  3 

2.2  The  p-k  Method.  ..  . . 4 

2.3  Viscous  Damping  Model  for  Coupled  Motions .  5 

3.0  RESULTS  AND  DISCUSSIONS . 6 

3.1  Aerodynamic  Coefficients .  6 

3.2  Flutter  Analysis  of  NACA64A006  Airfoil . 7 

3.3  Response  of  Airfoil  to  External  Excitations .  7 

4.0  CONCLUSIONS .  8 

5.0  REFERENCES . 9 

Table  1  Distances  from  Leading  Edge  for  Cp  Plotted  in  Figure  21 .  11 


ILLUSTRATIONS 


Figure  Page 

1  Two-Degree-of-Freedom  Airfoil  Motion .  13 

2  Variation  of  Lift  Coefficient  Due  to  Plunge  With  Reduced  Frequency  for  a 

NACA64A006  Airfoil  at  M  =  0.8 .  14 

3  Variation  of  Moment  Coefficient  Due  to  Plunge  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.8 .  15 

4  Variation  of  Lift  Coefficient  Due  to  Pitch  With  Reduced  Frequency  for  a 

N AC A64A006  Airfoil  at  M  =  0.8 . . .  16 

5  Variation  of  Moment  Coefficient  Due  to  Pitch  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.8 .  17 

6  Variation  of  Lift  Coefficient  Due  to  Plunge  With  Reduced  Frequency  for  a 

N AC A64A006  Airfoil  at  M  =  0.85 .  18 

7  Variation  of  Moment  Coefficient  Due  to  Plunge  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.85 . .  19 


(iv) 


ILLUSTRATIONS  (Cont’d) 


Figure  Page 

8  Variation  of  Lift  Coefficient  Due  to  Pitch  With  Reduced  Frequency  for  a 

NACA64A006  Airfoil  at  M  =  0.85 .  20 

9  Variation  of  Moment  Coefficient  Due  to  Pitch  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.85 .  21 

10  Variation  of  Lift  Coefficient  Due  to  Plunge  With  Reduced  Frequency  for  a 

NACA64A006  Airfoil  at  M  =  0.875  .  ...  22 

11  Variation  of  Moment  Coefficient  Due  to  Plunge  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.875  .  23 

12  Variation  of  Lift  Coefficient  Due  to  Pitch  With  Reduced  Fri  qutncy  for  a 

NACA64A006  Airfoil  at  M  =  0.875  .  24 

13  Variation  of  Moment  Coefficient  Due  to  Pitch  With  Reduced  Frequency 

for  a  NACA64A006  Airfoil  at  M  =  0.875  .  25 

14  Variation  of  Lift  Coefficient  Due  to  Plunge  With  Mach  number  for  a 

NACA64A006  Airfoil  at  kc  =  0.1  and  0.2  . .  26 

15  Variation  of  Moment  Coefficient  Due  to  Plunge  With  Mach  Number  for  a 

NACA64A006  Airfoil  at  kc  =  0.1  and  0.2  . . .  27 

16  Variation  of  Lift  Coefficient  Due  to  Pitch  With  Mach  Number  for  a 

NACA64A006  Airfoil  at  kc  =  0.1  and  0.2 . .  28 

17  Variation  of  Moment  Coefficient  Due  to  Pitch  With  Mach  Number  for  a 

NACA64A006  Airfoil  at  kc  =  0.1  and  0.2 .  29 

18  Cp  Distribution  for  Plunging  Motion  on  the  Upper  ( - )  and  Lower  ( - •) 

Surfaces  of  a  NACA64A006  Airfoil  at  M  =  0.85  and  kc  =  0.5 .  30 

19  Cp  Distrib  >n  for  Plunging  Motion  on  the  Upper  ( — ■ — )  and  Lower  ( - ) 

Surfaces  of  a  NACA64A006  Airfoil  at  M  =  0.85  and  kc  =  2 .  31 

20  Shock  Wave  Position  for  Plunging  Oscillation  of  a  NACA64A006  Airfoil  at 

M  =  0.85,  kc  =  0.5  and  2 .  32 

21  Cp  Variation  With  Time  on  Upper  Surface  of  a  NACA64A006  Airfoil  at 

M  =  0.85,  kc  =  0.5  and  2.  (The  distances  from  leading  edge  where  Cp 

are  computed  are  given  in  Table  1) .  33 

22  Variation  of  \J/bu>a,  co/cuft  and  kt.  With  \i  for  a  NACA64A006  Airfoil  at 

M  =  0.85 .  34 

23(a)  Comparison  of  Damping  Ratio  and  cj,/coq  With  U/buQ  Between  U-gand 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  n  =  50 .  35 


(v) 


Figure 

23(b) 


ILLUSTRATIONS  (Cont’d) 


Page 


■A'. 


a 

B 


Comparison  of  Damping  Ratio  and  co2/(ua  With  U/bcoa  Betveen  U-gand 

P-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  =  50 .  36 

24(a)  Comparison  of  Damping  Ratio  and  With  U/bcoa  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  n  =  15 .  37 

24(b)  Comparison  of  Damping  Ratio  and  o>2 /<Ja  With  U/bcoa  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  M  =  0.85  and  /i  =  75 .  38 

25(a)  Comparison  of  Damping  Ratio  and  o>,/a)a  With  U/bcoa  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  p  -  100  . .  39 


i 


25(b) 


Comparison  of  Damping  Ratio  and  cj2/ou a  With  U/bcoa  Between  U-g  and 
p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  /a  =  100  .... 


40 


I 


I 


26(a)  Comparison  of  Damping  Ratio  and  cji/u>a  With  U/bcaa  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  n  =  150 .  41 

26(b)  Comparison  of  Damping  Ratio  and  co2/coa  With  U/bcoa  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  p  =  150 .  42 

27(a)  Comparison  of  Damping  Ratio  and  a>[/tJa  With  U/ba>a  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  p  =  250  .  43 

27(b)  Comparison  of  Damping  Ratio  and  co2/coa  With  U/bu>a  Between  U-g  and 

p-k  Methods  for  a  NACA64A006  Airfoil  at  M  =  0.85  and  #1  =  250  .  44 

28  Variation  of  Amplitude  of  Displacement  Response  With  U/bcoa  and  to/coa 

for  an  Externally  Applied  Sinusoidal  Force .  45 

29  Variation  of  Phase  Angle  of  Displacement  Response  With  U/bcua  and  co/toa 

for  an  Externally  Applied  Sinusoidal  Force .  46 

30  Variation  of  Amplitude  of  Pitch  Response  With  U/bcoa  and  ce/cea  for  an 

Externally  Applied  Sinusoidal  Force .  47 

31  Variation  of  Phase  Angle  of  Pitch  Response  With  U/bwa  and  ce/cea  tor  an 

Externally  Applied  Sinusoidal  Force .  48 

32  Variation  of  Amplitude  of  Displacement  Response  With  U/bcoa  and  ce/c oa 

for  an  Externally  Applied  Sinusoidal  Moment  at  Elastic  Axis . .  49 

33  Variation  of  Phase  Angle  of  Displacement  Response  With  U/bcoa  and  ce/coa 

for  an  Externally  Applied  Sinusoidal  Moment  at  Elastic  A^is . .  50 

34  Variation  of  Amplitude  of  Pitch  Response  With  U/bcua  and  t  j/cua  for  an 

Externally  Applied  Sinusoidal  Moment  at  Elastic  Axis .  51 


-v- . 


v- 


(vi) 


ILLUSTRATIONS  (Cont’d) 


Figure  Page 

35  Variation  of  Phase  Angle  of  Pitch  Response  With  U/'bcoa  and  co/coa  for  an 

Externally  Applied  Sinusoidal  Moment  at  Elastic  Axis .  52 

36  Comparison  of  Amplitude  of  Displacemert  Response  Between  Viscous 

Damping  and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Force .  53 

37  Comparison  of  Phase  Angle  of  Displacement  Response  Between  Viscous 

Damping  and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Force .  54 

38  Comparison  of  Amplitude  of  Pitch  Response  Between  Viscous  Damping 

and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Force .  55 

39  Comparison  of  Phase  Angle  of  Pitch  Response  Between  Viscous  Damping 

and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Force .  56 

40  Comparison  of  Amplitude  of  Displacement  Response  Between  Viscous 

Damping  and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Moment 

at  Elastic  Axis .  57 

41  Comparison  of  t'base  Angle  of  Displacement  Response  Between  Viscous 

Damping  and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Moment 

at  Elastic  Axis .  58 

42  Comparison  cf  Amplitude  of  Pitch  Response  Between  Viscous  Damping 

and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Moment  at 

Elastic  Axis . .  .  59 

43  Comparison  of  Phase  Angle  of  Pitch  Response  Between  Viscous  Damping 

and  p-k  Methods  for  an  Externally  Applied  Sinusoidal  Moment  at 

Elastic  Axis .  60 


SYMBOLS 

Symbol  Definition 

ah  non-dimensional  distance  measured  from 

b  semi-chord  of  airfoil 

c  chord 

Cp  pressure  coefficient 

Clh  lift  coefficient  due  to  plunge 

Ct(>  lift  coefficient  due  to  pitch 

Cmh  moment  coefficient  due  to  plunge 


midchord  to  elaslic  axis 

i 


(vii) 


Symbol 

g 

h 

kc 

m 

M 

P 

P 

q 

Qh 


q; 

Q* 

q: 

t 

u 

xa 


6v, 


7 

‘-‘a 

a 

7 

r 

fh 

fa 

X 


SYMBOLS  (Cont’d) 

Definition 

moment  coefficient  due  to  pitch 
structural  damping  coefficient 
plunge  displacement 
reduced  frequency  coc/U 
mass  of  airfoil  per  unit  span 
free  stream  Mach  number 
complex  eigenvalue 
non-dimensional  complex  eigenvalue 
dynamic  pressure 
total  forces 

externally  applied  force 

total  moment 

externally  applied  moment 

radius  of  gyration  about  elastic  axis 

time 

free  stream  velocity 

non-dimensional  distance  measured  from  elastic  axis  to  center  of  mass 
defined  in  Equation  34 
defined  in  Equation  35 
pitch  angle 

damping  ratio  defined  in  Equation  18 
damping  ratio  defined  in  Equation  31 
viscous  damping  for  plunging  motion 
viscous  damping  for  pitching  motion 
flutter  eigenvalue  defined  in  Equation  9 


(viii) 


Symbol 

P 

* 

P 


u> 

<*>h 

“a 


SYMBOLS  (Cont’d) 
Definition 

airfoil-air  mass  ratio,  m/jrpb2 
non-dimensional  displacement  h/c 
air  density 
frequency 

uncoupled  bending  frequency 
uncoupled  torsional  frequency 


(ix) 


A  STUDY  OF  TRANSONIC  FLUTTER  OF  A  TV.O-DIMENSIONAL  AIRFOIL 


USING  THE  U-g  AND  p-k  METHODS 


1.0  INTRODUCTION 

Two  and  three-degree-of-freedom  transonic  flutter  of  two-dimensional  conventional  and 
supercritical  airfoils  has  been  the  subject  of  numerous  investigations  in  recent  years  (Refs.  1-12).  The 
drop  in  flutter  velocity  at  transonic  Mach  number  is  an  important  factor  in  design  of  modem  airfoils. 

In  the  past,  before  unsteady  transonic  computation  codas  were  available,  the  flutter  speed  in  the 
transonic  range  was  estimated  by  extrapolating  the  subsonic  and  supersonic  branches  of  the  flutter 
curves.  Among  the  methods  for  predicting  unsteady  transonic  flow  past  oscillating  airfoil  (Refs.  13-19), 
those  based  on  the  transonic  small  disturbance  equation  are  most  efficient,  and  hence  are  widely  used 
in  flutter  calculations  of  two-dimensional  airfoils. 

Two  approaches  Lave  been  employed  in  the  analysis  of  transonic  flutter.  The  first  utilizes 
the  conventional  U-g  method  (Ref.  20).  The  aerodynamic  lift  and  moment  coefficients  are  determined 
from  transonic  flow  calculations  and  the  structural  equations  are  soi.ed  for  the  flutter  velocity 
(Refs.  1-5).  Airfoils  that  have  been  studied  include  NACA64A006,  NACA64A010  ana  MBB  A-3  with 
variations  in  parameters  such  as  mass  ratio,  center  of  mass  location,  frequency  ratio,  etc. 

The  second  approach  determines  the  aeroelastic  responses  by  coupling  aerodynamic 
computation  codes  with  structural  equations  of  motion  using  a  variety  of  numerical  integration 
techniques  (Refs.  6-12).  Such  an  approach  was  first  taken  by  Ballhaus  and  Goorjian  (Ref.  6)  who 
carried  out  the  aeroelastic  response  of  a  NACA64A006  airfoil  with  a  single  degree-of-freedom  in 
pitch.  Extensions  of  this  procedure  to  two  and  three-degree-of-freedom  and  other  airfoils  such  as 
NACA64A010  and  MBB  A-3  have  been  reported  by  Yang  et  al.  (Refs.  8-10). 

The  original  finite-difference  scheme  for  solving  the  transonic  small  disturbance  equation, 
called  LTRAN2  by  Ballhaus  and  Goorjian  (Ref.  16),  is  accurate  only  at  low  reduced  frequencies 
(kc  <  0.075).  Houwink  and  van  der  Vooren  (Ref.  17)  improved  this  code  by  retaining  the  time- 
derivative  terms  in  the  boundary  and  auxiliary  conditions.  Their  modified  code  called  LTRAN2-NLR 
claimed  accuracy  up  to  kc  =  0.4.  Similar  modifications  have  also  been  carried  out  by  Couston  and 
Angelini  (Ref.  18)  at  ONERA.  Further  improvements  have  been  reported  by  Rizzetta  and  Chin 
(Ref.  19)  who  retained  the  second  time -derivative  term  of  the  velocity  potential  in  the  governing 
equation  and  presented  results  for  a  NACA64A10  airfoil  with  flap  oscillations  up  to  reduced 
frequency  of  5.  Isogai  (Ref.  5)  developed  a  similar  computation  code,  and  claimed  accuracy  up  to 
kc  =  0.25  for  the  entire  transonic  Mach  number  range  (from  subcritical  Mach  number  to  above  Mach 
number  1).  Similar  improvements  to  the  Couston’s  and  Angelini’s  code  (Ref.  18)  have  been  carried 
out  at  ONERA,  and  this  code  has  been  made  available  to  NAE.  Minor  modifications  (Ref.  21)  have 
been  made  and  these  included  addition  of  some  graphics  for  the  output  data.  The  use  of  these  various 
modified  or  improved  unsteady  t'ansonic  aerodynamic  codes  in  flutter  analysis  can  be  found  in 
References  1-13. 

The  U-g  method  gives  critical  flutter  speeds  in  agreement  with  other  methods  such  as  the 
traditional  British  approach  with  lined-up  frequency  parameters,  but  overestimates  the  relative  damp¬ 
ing  ratio  at  other  speeds  (Ref.  22).  The  use  of  Frueh’s  and  Miller’s  (Ref.  23)  correction  formula  for 
damping  gives  useful  values  of  the  decay  factors,  but  Woodcock  and  Lawrence  (Ref.  24)  found  the 
U-g  method  is  still  inferior  to  the  British  method.  Hassig  (Ref.  25)  developed  a  method  which  gives 
damping  values  in  excellent  agreement  with  the  British  method.  The  p-k  method,  as  called  by  Hassig 
(Ref.  25),  iterates  for  the  zeros  of  the  flutter  determinant  using  values  of  the  aerodynamics  interpo¬ 
lated  from  given  values  at  a  set  of  frequency  parameters.  This  method  assumes  that  for  sinusoidal 
motions  with  slowly  varying  amplitude,  the  aerodynamic  forces  can  be  approximated  by  those  based 
on  constant  amplitude. 


The  use  of  time  marching  techniques  which  couple  the  aerodynamics  with  the  structural 
equations  of  motion  should  g;ve  accurate  resulis  for  the  critical  flutter  speeds  and  damping  ratios  at 
other  speeds.  However,  this  method  is  time  consuming  and  is  not  suited  for  flutter  analysis  where 
many  parameters  such  as  mass  ratio,  frequency  ratio,  etc.  are  varied. 

The  first  part  of  this  report  presents  aerodynamics  data  of  a  NACA64A006  airfoil  oscillating 
in  p  h  and  plunge  using  the  high  frequency  version  of  ON  ERA  unsteady  transonic  code.  The  results 
are  compared  with  those  obtained  by  Yang  and  Cher.  (Ref,  9)  employing  LTRAN2-NLR.  To  the  best 
of  the  author’s  knowledge,  data  from  the  ONERA  code  have  not  been  reported  by  ONERA  in  the 
open  literature.  They  are,  therefore,  presented  in  somewhat  detailed  manner  in  this  report.  From  the 
aerodynamics  obtained  from  the  ONERA  code,  the  U-g  method  is  used  to  compute  critical  flutter 
speeds  and  damping  ratios  for  a  two-degree-of-freedom  NACA64A006  airfoil  using  Frueh’s  and 
Miller’s  formula  (Ref.  23).  The  results  using  the  p-k  method  are  compared  with  those  from  the  U-g 
method.  Response  curves  from  the  p-k  method  and  those  using  a  viscous  damping  model  for  coupled 
motions  are  also  included  for  forced  sinusoidal  motion. 


2.0  FLUTTER  ANALYSIS  OF  TWO-DEGREE-OF-FREEDOM  AIR  FOIL  MOTION 

Figure  1  shows  the  notations  used  in  the  analysis  of  a  two-degree-of-freedom  motion  of  an 
airfoil  oscillating  in  pitch  and  in  plunge.  The  bending  deflection  is  denoted  by  h,  positive  in  the 
downward  direction,  a  is  the  pitch  angle  about  the  elastic  axis,  positive  with  the  nose  up.  The  elastic 
axis  is  located  at  a  distance  ahb  from  the  midchord,  while  the  mass  center  is  located  at  a  distance  x0b 
from  the  elastic  axis.  Both  distances  are  positive  when  measured  towards  the  trailing  edge  of  the 
axrfoil.  The  aeroelastic  equations  of  -notion  have  been  derived  by  Fung  (Ref.  20)  and  can  bt  written 
as: 


2  t  = 


£  +  xaor  +  wj|  $ 


Qh 

mb 


Qa 


mbz 


(1) 


(2) 


where  £  =  h/b  is  the  non-dimensional  displacement,  m  is  the  mass  per  unit  span  of  ,ne  airfoil,  coh,  coa 
are  the  uncoupled  plunging  and  pitching  frequency  respectively,  r0  is  the  radius  of  gyration  about 
the  elastic  axis,  and  Qh  and  Qa  are  the  total  forces  r.nd  moments  acting  about  the  elastic  axis 
respectively.  In  the  absence  of  externally  applied  forces  and  moments,  Qh  and  Qa  can  be  written  as: 


Qh  =  -qc  (cCh  |  +  C5q 

Qa  =  qc2(cmh|  +  Cmao) 


(3) 


(4) 


where  q  is  the  dynamic  pressure,  Cc  and  Cm  are  the  lift  and  moment  coefficients  with  the  subscripts 
h  and  a  denoting  unit  plunging  and  pitching  motionr  respectively. 


2.1  The  U-g  Method 


In  the  U-g  method,  a  structural  damping  coefficient  g  is  introduced  into  Equations  (1) 
and  (2)  by  multiplying  the  third  term  of  the  two  equations  by  the  factor  (1  +  ig).  For  harmonic 
oscillations,  £  and  a  can  be  written  in  the  following  form: 


*  =  *0  eitJ* 

a  =  a0  eilJt 

Equations  (1)  and  (2)  can  be  expressed  as: 


where  p  =  m/rrb  2p  is  the  airfoil-air  mass  ratio,  kc  =  ojc/U  is  the  reduced  frequency  and 
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The  complex  eigenvalue  X  can  be  solved  from  the  following  quadratic  equation 
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Hassig  (Ref.  25),  quoting  from  Frueh  and  Miller  (Ref.  23),  gives  an  expression  for  the 
damping  ratio  which  is  written  in  the  present  notations  as: 
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2.2  The  p-k  Method 

Equations  (5)  and  (6)  can  be  written  in  the  following  form: 


t  =  *p' 


(19) 


J.-.- 

•'S 


a  =  a0  ept 


where  p  =  0  +  ico,  and  0  is  the  damping  factor.  Substituting  into  Equations  (1)  to  (4)  gives: 
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where  p  =  —  •  p,  and  can  be  solved  from  the  following: 


Dp4  +  Ep2  +  F  =  0  (23; 

where  D  =  r2  -  x2  (24 

E  =  d  +  ar2  -  xa  (c  +  b)  (25 

F  =  ad  -  be  (26 


The  damping  ratio  f  is  given  as: 


f  =  - 


V/32  +  to2 


(31 


2.3  Viscous  Damping  Model  for  Coupled  Motions 


In  the  study  of  tbe  dynamic  response  of  a  v-lug  to  random  loads  or  buffeting,  a  viscous 
damping  force  is  often  used  in  the  structural  equations  of  motion  in  place  of  the  aerodynamic 
damping  force,  and  the  damping  ratio  is  obtained  from  flutter  calculations  (see,  for  example,  Ref.  26). 
Introducing  a  viscous  damping  term  in  Equations  (1)  and  (2).  and  assuming  harmonic  motion  of  the 
form  given  in  Equations  (5)  and  (6),  the  structural  equations  of  motion  are: 
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v/here 
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(35) 

fh  and  fa  are  the  damping  ratios  for  plunging  and  pitching  motions.  Equations  (32)  and  (33)  can  be 
used  to  study  the  response  of  the  airfoil  to  external  loads. 


3.0  RESULTS  AND  DISCUSSIONS 
3.i  Aerodynamic  Coefficients 

In  this  study  only  two-degree-of -freedom  motions  in  pitch  and  plunge  are  considered.  The 
aerodynamic  coefficients  presented  are  those  for  lift  and  moment  due  to  unit  plunge  displacement 
and  pitch  rotation  respectively,  namely,  C„h,  Cmh,  CCo,  and  Cma.  The  computations  are  performed 
using  the  ONERA  improved  version  of  Couston’s  and  Angelinas  (Ref.  18)  unsteady  transonic  aero¬ 
dynamics  code.  The  lesults  for  a  NACA64A006  airfoil  are  presented  herein  in  a  somewhat  detailed 
manner  since  it  does  not  appear  that  computations  using  this  computer  code  have  been  compared 
with  other  transonic  codes  and  reported  in  the  open  literature. 

In  Figures  2  to  13,  CBh,  Cmh,  CB(>,  and  Cma  are  plotted  versus  kc  for  Mach  numbers 
M  =  0.8,  0.85  and  0.875.  The  amplitudes  for  the  plunging  and  pitching  motions  are  h/c  =  0.04  and 
a  =  0.1  deg.  respectively,  and  the  pitch  axis  is  located  at  the  1/4  chord  point.  The  results  are  computed 
up  to  a  value  of  kc  =2.  Shown  also  in  the  figure.,  are  computations  by  Yang  and  Chen  (Ref.  9)  using 
LTRAN2-NLR  which  is  claimed  to  be  accurate  up  to  kc  =  0.4.  The  results  from  the  ONERA  code 
agree  quite  well  with  LTRAN2-NLR  up  to  and  in  some  cases  exceeding  kc  =  0.4,  except  for  the  Cm 
curves  (especially  the  imaginary  part)  at  M  =  0.85.  At  higher  frequencies,  LTRAN2-NLR  is  not  “ 
accurate  since,  unlike  the  ONERA  code,  the  second  time  derivative  of  the  velocity  potential  term  has 
not  been  included  in  the  governing  equation. 

The  four  aerodynamic  coefficients  are  plotted  versus  Mach  number  in  Figures  14  to  17  for 
reduced  frequency  kc  =  0.1  and  0.2.  It  is  seen  that  at  M  =  0.85  there  are  large  changes  in  the  aero¬ 
dynamic  coefficients,  especially  the  moment  coefficients.  This  Mach  number  coincides  with  the 
transonic  dip  observed  for  the  NACA64A006  airfoil  by  Yang  et.  al.  (Refs.  1,  2).  At  M  =  0.7  and  0.8 
the  flow  past  the  airfoil  contains  no  shock  waves.  At  M  =  0.85,  Fig.  18  shows  that  a  shock  wave  is 
present  for  part  of  a  plunge  cycle  at  kc  =  0.5.  The  results  are  typical  up  to  a  value  of  kc  approximately 
1.4,  above  which  a  shock  wave  is  present  for  the  full  cycle  as  shown  in  Figure  19  for  kc  =  2.  The 
shock  positions*  for  the  third  and  fourth  cycles  of  computations  are  shown  in  Figure  20  for  kc  =0.5 
and  2.0.  It  is  seen  that  at  low  frequencies,  the  shock  movements  are  much  larger  than  those  at  higher 
frequencies. 

For  plunging  motions,  Cp  plots  for  the  upper  surface  of  the  airfoil  at  M  =  0.85  for  the  third 
and  fourth  cycles  of  computations  are  shown  in  Figure  21  for  kc  =0.5  and  2.0.  The  numbers  from  3 
to  29  marked  on  the  curves  represent  computation  grid  points,  and  their  distances  from  the  leading 
edge  non-dimensionalized  wj.t.  the  chord  are  given  in  Table  1.  Because  of  the  interrupted  shock  wave 


'The  shock  positions  are  determined,  according  to  Reference  21,  from  grid  points  where  ACp/Ax  is  a  maximum,  while  also  satisfying 
the  criterion  ACp/Ax  >  2 


motions  at  low  frequt  ncies,  shown  for  example  in  Figure  21  (kc  =  0.5),  the  pressures  in  the  region 
behind  the  shock  shov  irregular  oscillations.  At  higher  frequencies,  when  shock  waves  are  presented 
at  all  times  in  a  plunge  cycle  (Fig.  21,  kc  -  2.0),  the  pressures  behind  the  shock  oscillates  sinusoidally. 


Results  for  pitch  oscillations  similar  to  those  presented  in  Figures  18  to  21  are  not  reported 
in  this  report.  For  a  NACA64A010  airfoil,  Isogai  (Ref.  5)  gives  some  quite  detailed  results  for  pitching 
motion. 


3.2  Flutter  Analysis  of  NACA64A006  Airfoil 

Flutter  calculations  using  the  ONERA  unsteady  transonic  aerodynamics  code  are  compared 
in  Figure  22  with  those  obtained  by  Yang  et  al.  (Ref.  1)  using  an  implicit  scheme  (Ref.  6)  and  a 
relaxation  approach  (Refs.  14, 15)  for  the  aerodynamic  coefficients.  In  this  figure,  the  U-g  method  is 
used  and  the  flutter  speed  U/bcoa,  reduced  frequency  kc  and  frequency  ratio  cj/coa  are  plotted  against 
airfoil-air  mass  ratio  p.  The  Mach  number  is  0.85,  and  the  following  parameters  are  used:  xa  =  0.25, 
rQ  =  0.5,  wh  ju>a  =  0.2,  pitch  axis  =  0.25  (ah  =  -0.5),  which  are  the  ones  used  by  Yang  et  al.  (Ref.  1). 
This  figure  shows  lower  flutter  speeds  using  ONERA  code  than  those  from  Reference  6,  except  at  a 
mass  ratio  of  50.  The  present  results  can  be  taken  to  be  more  accurate  than  Reference  1  since  the 
aerodynamics  are  improved  by  taking  the  second  time  derivative  of  the  velocity  potential  in  the 
governing  equation.  The  aerodynamic  coefficients  are  in  good  agreement  with  LTRAN2-NLR  (see 
Figs.  6  to  9)  in  the  range  of  values  for  kc  between  0.1  and  0.2  where  flutter  speeds  are  calculated. 
Comparisons  of  the  aerodynamics  from  the  indicial  and  relaxation  methods  given  in  Reference  2  with 
either  the  LTRAN2-NLR  or  the  present  results  show  tne  relaxation  approach  yields  closer  agreement 
than  the  indicial  method  which  gives  much  poorer  results  especially  for  the  CCa  and  CWa  coefficients. 

The  effect  of  p  on  the  damping  and  frequency  ratios  are  shown  in  Figures  23  to  27  for 
M  =  0.85.  The  results  for  f  and  y  of  the  two  modes  are  computed  using  Equations  (18)  and  (31).  The 
critical  flutter  speeds  and  frequencies  obtained  from  both  methods  are  identical.  For  the  bending 
mode,  the  differences  in  subcritical  damping  and  frequency  ratios  between  the  p-k  and  U-g  methods 
decrease  for  increasing  p.  At  p  =  250,  these  two  methods  give  almost  identical  results.  For  the  torsion 
mode,  the  differences  between  f  and  y  increase  with  ju  while  the  differences  in  co2  /<uQ  decrease.  The 
fairly  good  agreement  between  the  two  methods  may  be  a  special  case  for  a  two-degree-of-freedom 
two-dimensional  airfoil  motion.  Results  for  subcritical  damping  from  Reference  24  show  poor  com¬ 
parison  between  the  U-g  and  the  British  method  with  lined-up  frequency  parameter,  which  has  in  turn 
been  demonstrated  by  Hassig  (Ref.  25)  to  be  in  good  agreement  with  the  p-k  method. 


3.3  Response  of  Airfoil  to  External  Excitations 


To  investigate  the  response  of  the  airfoil  as  the  flutter  speed  is  approached,  £  and  a  are 
solved  from  Equations  (1)  and  (2)  for  two  cases:  Qa  =  0  and  =  0,  where  Qa  and  are  the 
externally  applied  sinusoidal  excitation  moment  and  force  at  the  airfoil’s  elastic  axis  respectively. 
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in  Figures  28  to  35.  For  a  purely  translational  excitation  force,  the  displacement  amplitude  (Fig.  28) 
decreases  with  increasing  values  of  U/ba>Q  until  a  minimum  is  reached  around  U/bco0  =  2.5  to  3.0. 
Further  increase  in  U/bco0  results  in  a  rapid  growth  of  the  amplitude.  The  pitch  amplitude  (Fig.  30), 
shows  a  similar  behaviour  for  a>/cja  in  the  vicinity  of  the  uncoupled  bending  frequency.  The 
amplitude  near  the  uncoupled  torsional  frequency  decays  as  the  flutter  speed  is  approached. 
Figures  29  and  31  show  the  corresponding  phase  angles  for  displacement  and  pitch  rotation 
respectively. 


When  a  moment  is  applied  at  the  elastic  axis,  the  displacement  increases  rather  gradually 
initially  until  U/bcoQ  approaches  tne  flutter  speed  (Fig.  32).  Similar  behaviour  of  the  pitch  amplitude 
near  the  uncoupled  bending  frequency  is  shown  in  Figure  34.  The  peak  response  near  co/oJa  =  1 
decays  from  large  values  at  small  U/btoa  to  very  small  values  as  U/bcoa  approaches  the  critical  value. 
The  phase  angles  for  the  displacement  and  pitch  response  are  shown  in  Figures  33  and  35. 

In  Section  2.3  a  viscous  damping  model  far  the  airfoil  motion  is  given.  In  the  study  of 
dynamic  response  of  wing  to  external  loads,  such  as  buffeting,  often  the  aerodynamic  damping  is 
combined  with  the  structural  damping  to  yield  a  total  damping  proportional  to  the  velocity  (see,  for 
example,  Ref.  27).  The  damping  values  can  be  determined  either  experimentally  or  from  flutter 
calculations,  such  as  the  p-k  method.  Figures  36  to  43  show  the  displacement  and  pitch  response  of 
a  NACA64A006  airfoil  at  M  =  0.85  and  p  =  100.  The  values  of  fh,  and  coh /coQ  are  obtained  for 
U/bcoa  =  4  using  the  p-k  method.  Cu^/es  for  other  values  of  U/bcoa  are  quite  similar  to  those  shown 
in  these  figures.  For  a  translational  excitation  force  (Q'a  -  0),  Figures  36  to  39  show  the  viscous 
damping  and  p-k  methods  give  quite  similar  response  curves.  However,  when  the  external  excitation 
is  a  moment  applied  at  the  elastic  axis  (Q^  =  0),  Figures  40  and  42  show  the  large  response  obtained 
near  the  uncoupled  bending  frequency  from  the  p-k  method  cannot  be  duplicated  using  the  viscous 
damping  model,  even  though  the  response  for  pitching  motion  near  the  uncoupled  torsional  frequency 
agrees  fairly  well.  The  phase  curves  for  the  displacement  show  a  difference  of  nearly  180  degrees 
between  the  two  methods,  while  those  for  the  pitch  rotation  are  reasonably  good  except  near  the 
uncoupled  bending  frequency. 


4.0  CONCLUSIONS 

The  improved  version  of  ONERA  unsteady  transonic  aerodynamics  code,  which  include 
the  second  time  derivative  term  of  the  velocity  potential  in  the  governing  equation,  has  been  used  to 
compute  aerodynamic  coefficients  for  plunging  and  pitching  motions  of  a  NACA84A006  airfoil  at 
various  Mach  numbers.  The  results  agree  quite  well  with  LTRAN2-NLR  up  to  and  in  some  cases 
exceeding  kc  =  0.4,  exce  >t  for  the  Cma  curves  (especially  the  imaginary  part)  at  M  =  0.85.  This 
Mach  number  coincides  with  the  transonic  dip  and  large  changes  in  the  aerodynamic  coefficients, 
especially  the  moment  coefficients,  are  observed.  Results  of  shock  motions  are  given  for  the  airfoil 
oscillating  in  plunge.  At  M  =  0.85  and  kc  below  1.4,  computations  show  that  a  shock  wave  is  present 
only  for  part  of  a  plunge  cycle.  The  shock  motions  are  fairly  large  and  the  pressures  in  the  region 
behind  the  shock  show  irregular  oscillations. 

Flutter  speeds  for  a  NACA64A006  airfoil  using  the  U-g  method  and  aerodynamic 
coefficients  from  ONERA  code  are  lower  than  those  reported  by  Yang  et  al.  using  LTRAN2.  The 
comparisons  are  made  at  M  =  0.85  for  various  airfoil -air  mass  ratios.  The  p-k  method  gives  flutter 
speeds  and  kc  values  identical  to  those  from  the  U-g  method.  The  computed  subcritical  damping 
ratios  at  M  =  0.85  for  different  n  show  the  U-g  method  using  Frueh’s  and  Miller’s  formula  gives 
results  quite  close  to  the  p-k  method,  especially  for  large  values  of  pt.  The  comparison  is  always  very 
close  for  small  values  of  U/bcoa  irrespect  of  pi.  The  fairly  good  agreement  between  these  two  methods 
may  be  a  special  case  for  a  two-degree-of-freedom  two-dimensional  airfoil  performing  plunging  and 
pitching  oscillations.  Results  for  more  complicated  cases  show  that  subcritical  damping  ratios  obtained 
from  the  U-g  method  agree  poorly  with  the  British  method  with  lined-up  frequency  parameters,  which 
has  in  turn  been  demonstrated  to  be  in  good  agreement  with  the  p-k  method. 

The  p-k  method  can  be  used  to  study  the  response  of  an  airfoil  to  externally  applied 
excitation  forces  and  moments.  Two  examples  are  given  for  either  a  translational  force  or  a  moment 
ap»°ied  at  the  elastic  axis.  Using  the  computed  values  of  damping  ratios  fh  and  for  the  bending 
a  id  torsion  modes  together  with  the  ratio  oj/coa,  the  response  using  a  viscous  damping  model  has 
been  determined  and  compared  with  the  p-k  method.  It  is  shown  that  for  translational  excitation 
force,  the  response  curves  for  the  two  methods  are  in  good  agreement.  However,  when  the  external 
excitation  is  a  moment  applied  at  the  elastic  axis,  the  large  response  obtained  near  the  uncoupled 
bending  frequency  from  the  p-k  method  cannot  be  duplicated  using  the  viscous  damping  model,  even 
though  the  response  for  pitching  motion  near  the  uncoupled  torsional  frequercy  agrees  fairly  well. 
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TABLE  1 

DISTANCES  FROM  LEADING  EDGE  FOR  Cp  PLOTTED  IN  FIGURE  21 


Position 

x/c 

Position 

x/c 
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0.04029 

17 

0.5769 
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0.6712 
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0.1410 
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0.7546 
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0.2079 
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0.8269 
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0.2857 
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0.8883 
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0.3745 
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0.9386 
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0.4744 
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FIG.  1:  TWODEGREE-OF-FREEDOM  AIRFOIL  MOTION 
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FIG.  2:  VARIATION  OF  LIFT  COEFFICIENT  DUE  TO  PLUNGE  WITH  REDUCED 
FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.8 


FIG.  5:  VARIATION  OF  MOMENT  COEFFICIENT  DUE  TO  PITCH  WITH  REDUCEO 
FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.3 


FIG.  6:  VARIATION  OF  LIFT  COEFFICIENT  DUE  TO  PLUNGE  WITH  REDUCED 
FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85 
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FIG.  7:  VARIATION  OF  MOMENT  COEFFICIENT  DUE  TO  PLUNGE  WITH  REDUCED 
FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85 
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FIG.  8:  VARIATION  OF  LIFT  COEFFICIENT  DUE  TO  PITCH  WITH  REDUCED 
FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85 
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FIG.  13:  VARIATION  OF  MOMENT  COEFFICIENT  DUE  TO  PITCH  WITH  REDUCED 


FREQUENCY  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.875 
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FIG.  16:  VARIATION  OF  LIFT  COEFFICIENT  DUE  TO  PITCH  WITH  MACH  NUMBER 
FOR  A  NACA64A006  AIRFOIL  AT  kc  =  0.1  AND  0.2 
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FIG.  17:  VARIATION  OF  MOMENT  COEFFICIENT  DUE  TO  PITCH  WITH  MACH 
NUMBER  FOR  A  NACA64A006  AIRFOIL  AT  k  =0.1  AND  0.2 
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FIG.  20:  SHOCK  WAVE  POSITION  FOR  PLUNGING  OSCILLATION  OF  A  NACA64A006 

AIRFOIL  AT  M  =  0.85,  kc  =  0.5  AND  2 


FIG.  22:  VARIATION  OF  U/bwo,  oj/wa  AND  kc  WITH  ju  FOR  A  NACA64A006 

AIRFOIL  ATM  =  0.85 


FIG.  23(a):  COMPARISON  OF  DAMPING  RATIO  AND  <j  1/wo  WITH  U/bcoQ  BETWEEN 
U-g  AND  p-k  METHODS  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85  AND  n  =  50 


FIG.  23(b):  COMPARISON  OF  DAMPING  RATIO  AND  u2,ua  W,TH  U/btJ*  BETWEEN 
U-g  AND  p-k  METHODS  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85  AND  n  =  50 
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FIG.  25(b):  COMPARISON  OF  DAMPING  RATIO  AND  <*>2h\  WITH  U/bw,  BETWEEN 
U-g  AND  p-k  METHODS  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85  AND  p  =  100 
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FIG.  26(b):  COMPARISON  OF  DAMPING  RATIO  AND  u2/ua  WITH  U/bco,  BETWEEN 
U-g  AND  p-k  METHODS  FOR  A  NACA64A006  AIRFOIL  AT  M  =  0.85  AND  n  =  150 


o 


FIG.  34:  VARIATION  OF  AMPLITUDE  OF  PIT^H  RESPONSE  WITH  U/bcoa  AND  u/ua 
FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL  MOMENT  AT  ELAStIc  AXIS 
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FIG.  35:  VARIATION  OF  PHASE  ANGLE  OF  PITCH  RESPONSE  WITH  U/b<Ja  AND 
oj/cj  FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL  MOMENT  AT  ELASTIC  AXIS 
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FIG.  36:  COMPARISON  OF  AMPLITUDE  OF  DISPLACEMENT  RESPONSE  BETWEEN  VISCOUS 
DAMPING  AND  p-k  METHODS  FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL  FORCE 
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FIG.  39:  COMPARISON  OF  PHASE  ANGLE  OF  PITCH  RESPONSE  BETWEEN  VISCOUS 
DAMPING  AND  p-k  METHODS  FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL  FORCE 


FIG.  40:  COMPARISON  OF  AMPLITUDE  OF  DISPLACEMENT  RESPONSE  BETWEEN  VISCOUS 
DAMPING  AND  p-k  METHODS  FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL  MOMENT  AT 

ELASTIC  AXIS 


0.00  0.2S 


0.50 


0.75 


1.00 


1.25  1.50  1.75  2. CO 


FIG.  42:  COMPARISON  OF  AMPLITUDE  OF  PITCH  RESPONSE  BETWEEN  VISCOUS 
DAMPING  AND  p-k  METHODS  FOR  AN  EXTERNALLY  APPLIED  SINUSOIDAL 

MOMENT  AT  ELASTIC  AXIS 
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SUMMARY/SOMMAIRE 

Transonic  flutter  of  a  NACA64A006  airfoil  undergoing  plunging 
and  pitching  oscillations  is  studied  using  the  U-g  and  p-k  methods.  The 
aerodynamic  coefficients  are  calculated  using  an  improved  version  of  an 
ONERA  unsteady  transonic  aerodynamics  code  which  include  the  second 
time  derivative  term  of  the  velocity  potential  in  the  governing  equation. 
Comparisons  with  LTRAN2-NLR  show  good  agreement  up  to  and  in  some 
cases  exceeding  kc  =  0.4,  except  for  the  pitching  moment  curves  at  the 
transonic  dip  Mach  number  of  0.85.  All  flutter  results  are  presented  for 
M  =  0.85.  The  p-k  method  gives  flutter  speeds  identical  to  those  from  the 
U-g  method.  Subcritical  damping  ratios  using  the  U-g  method  with  Frueh’s 
and  Miller’s  damping  formula  are  quite  close  to  those  obtained  from  the  p-k 
method,  especially  for  large  values  of  the  airfoil-air  mass  ratio.  Response  of 
the  airfoil  to  externally  applied  forces  and  moments  is  studied  using  the  p-k 
method  and  a  viscous  damping  model  for  coupled  motions. 


